High fidelity quantum memory via dynamical decoupling: theory and experiment 



Xinhua Peng 1 , Dieter Suter 1 , Daniel A. Lidar 2 

1 Department of Physics, Technische Universitdt Dortmund, Germany 
2 Departments of Electrical Engineering, Chemistry, and Physics, 

Center for Quantum Information Science & Technology, 
University of Southern California, Los Angeles, CA 90089, USA 

Quantum information processing requires overcoming decoherence - the loss of "quantumness" due to the 
inevitable interaction between the quantum system and its environment. One approach towards a solution is 
quantum dynamical decoupling - a method employing strong and frequent pulses applied to the qubits. Here we 
report on the first experimental test of the concatenated dynamical decoupling (CDD) scheme, which invokes 
recursively constructed pulse sequences. Using nuclear magnetic resonance, we demonstrate a near order of 
magnitude improvement in the decay time of stored quantum states. In conjunction with recent results on high 
fidelity quantum gates using CDD, our results suggest that quantum dynamical decoupling should be used as 
a first layer of defense against decoherence in quantum information processing implementations, and can be a 
stand-alone solution in the right parameter regime. 



I. INTRODUCTION 



One of the main difficulties in the implementation of quan- 
tum information processing is the susceptibility of quantum 
systems to interactions with their environment. Any such 
uncontrolled interaction degrades the quantum information 
stored in the system; in the present context, this process is 
known as decoherence. 1 When a quantum system performs 
an information processing task, decoherence causes computa- 
tional errors, which leads to the eventual loss of any quan- 
tum advantage in information processing. 2 Dynamical de- 
coupling (DD) is a form of quantum error suppression that 
modifies the system-environment interaction so that its over- 
all effects are very nearly self-canceling, thereby decoupling 
the system evolution from that of the decoherence-inducing 
environment. 3 4 DD has been primarily studied as a special- 
ized quantum control technique for state preservation, or 
"quantum memory". In this setting it has been demonstrated 
theoretically that DD is capable of generating long-lived and 
robust quantum states. Three recently introduced techniques 
have figured most prominently in this effort: "randomized 
DD" (RDD), 5 "concatenated DD" (CDD), 6 and "Uhrig DD" 
(UDD). 7 RDD, which involves randomly selected pulse types 
applied at regular intervals, works best for strongly time- 
dependent environments and in the long time limit. 8 9 UDD, 
which involves an optimization of pulse intervals for a fixed 
pulse type, is provably optimal for the restricted case of a sin- 
gle qubit with diagonal coupling to the environment. 10 11 UDD 
has been the subject of extensive recent experimental tests. 
Following some notable successesP^ the general picture 
emerging is that UDD performance depends sensitively on 
the existence of a sharp high-frequency cutoff; in the absence 
of such a cutoff simpler DD schemes tend to result in better 
performance PK] UDD was recently extended to "quadratic 
DD" (QDD) which applies to general environments, but is 
still limited to only a single qubitPS Very recently a gener- 
alization of QDD to multiple qubits was proposed, which in- 
volves a higher order nesting of UDD sequences!^ CDD is 
a method for constructing pulse sequences with a recursive 
structure. While CDD requires longer pulse sequences than 



UDD and QDD, there is strong theoretical evidence that it can 
be successfully combined with quantum computation and in- 
corporated into quantum logic circuits p322 p or this reason we 
focus here on CDD, and present the first experimental evi- 
dence, using nuclear magnetic resonance (NMR) techniques, 
for the theoretically predicted advantag e * 6 ! 9 * 20 " ^ of using CDD 
for quantum memory preservation over its common counter- 
part "periodic dynamical decoupling" (PDD). a We also pro- 
vide a theoretical analysis in support of these experimental 
results. Though our experimental results illustrate the effec- 
tiveness of CDD specifically for NMR, the CDD framework is 
generally applicable to any decohering quantum system. Our 
results lend support to the expectation that quantum dynam- 
ical decoupling will prove an indispensable tool in scalable 
quantum information processing. 



II. DYNAMICAL DECOUPLING AND CDD 

We consider a system described by a Hamiltonian 

H = Hs + H B + H SB , (1) 

where Hs acts on the system degrees of freedom, Hb on the 
environment and Hs b couples the system to the environment. 
To suppress errors, DD allows the system to evolve for some 
time before applying a control pulse to redirect or refocus 
the evolution toward the error- free ideal, continually repeat- 
ing this process until some total evolution has completed: 

DD[t/(r)] = P N U(r ) • • • P 2 U(r )P 1 U(r ) = U(Nt ), 

(2) 

where U(tq) = Uq(tq)B(to) represents the unitary evolu- 
tion of the combined system and environment, or bath, for a 
duration of length tq, decomposed so that Uq(tq) determines 
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the ideal, desired error-free evolution generated by a piece- 
wise constant system Hamiltonian H$, and B (to) is a uni- 
tary error operator acting jointly on the system and bath. We 
use the symbol U to denote evolution in the presence of DD 
pulses. For now, we implicitly assume that the pulses Pj are 
sufficiently fast as to not contribute to the total time of the 
evolution. This assumption is not essential, but is useful in 
simplifying our present discussion. The simplest example is 
quantum memory, where Uo(to) = Is is the identity opera- 
tion or "free-evolution" of the system, and B (to) represents 
the deviation from the ideal dynamics caused by the presence 
of a bath. In this case, our goal is to choose pulses so that 
DD[[/(r)] = Is <8> B, the identity acting on the system and 
an arbitrary pure bath operator B. DD schemes differ in pre- 
cisely how these pulses Pj are chosen, with the only common 
constraint that the following basic "decoupling condition" is 
met™ 

^P£H SB Pa=0. (3) 

OL 

To be concrete, we will suppose that the pulses P a G 
{I, X, Y, Z} are Pauli operators, though again this assump- 
tion only serves to keep the explanation simple and is not 
strictly necessary. CDD generates pulse sequences by recur- 
sively building on a base sequence (or any 
other ordered pair of Pauli operators from the set {X, Y, Z}), 
as follows. The sequence is initialized as 

CDD [tf(To)] = U(r ) = U (t )B(t ) = U (t ), (4) 
and higher levels are generated via the rule 

c7 n+1 (T n+1 ) = CDD n+1 [[/(T )] = 
Z[^ n (r n )]X[^ n (r n )]Z[^ n (r n )]X[£r n (r n )], (5) 

where r n = 4 n To. Thus, at each level the total duration of the 
CDD pulse sequence grows by a factor of 4, while the pulse 
interval remains fixed at its initial value To. Note that we are 
allowing for the possibility of some non-trivial information 
processing operation Uq{tq), as implemented in Ref.l2Tl The 
choice of the base sequence is motivated by the observation 
that it satisfies the "decoupling condition" ^ in the quantum 
memory setting Uo(to) = Is, under the dominant "1 -local" 
system-bath coupling term = Y,a=x, y ,z Y^j a j ® Bf, 
where crj = X, = Y, and crj = Z denote the Pauli ma- 
trices acting on system qubit j, and {B^} are arbitrary bath 
operators. (The next order "2-local" coupling would have 
terms such as cr^crf ® , etc.) For this reason the base 
sequence is sometimes called the "universal decoupler". Sim- 
ilarly, the most common pulse sequence used thus far in DD 
experiments is "periodic DD" (PDD), which generates pulse 
sequences by periodically repeating the universal decoupler 
base sequence Z[]X[]Z[}X[}\ 

PDB k [U(r )} = (PDD X [C/(T )]) fc = U k (4kr ), (6) 

where PDDi[C/(r )] = CDDi[C/(t )]. Note that in Eq. {6} 
Uk refers to the /c-fold repetition of single cycles of the se- 
quence, while in Eq. ^ U n describes the n-fold concatena- 
tion of identical cycles. 
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FIG. 1 : Experimental comparison of CDD vs PDD vs free evolution 
on a 300MHz NMR spectrometer. Shown is the signal decay of the 
13 C-spin coherence of Adamantane (C10H16, depicted in the inset) 
vs. total evolution time. The spin echo is measured by first creating 
the equal superposition state (|0) + |l))/\/2, then subjecting it to 
either PDD (red squares) or CDD (blue circles) pulse sequences, and 
finally measuring the free induction decay. The first point shows the 
amplitude of the free induction decay (FID) signal, measured after a 
pulse interval of to = 15/xs [i.e., CDDo, Eq. j4|]. The second point 
corresponds to CDDi = PDD. The third and the fourth points 
correspond to CDD 2 and CDD 3 , respectively. The width of the 
7r-pulses was 8 = 10.52/xs. For reference we also show the free 
evolution, without any pulses (triangles). 



The use of periodic pulse sequences has a long history, 
starting with some of the earliest spin-echo experiments 2 ^. In- 
deed, it is natural to try to design a "good" pulse sequence and 
then repeat it over and over again so as to cover a desired total 
experiment time T. CDD breaks with this intuition and tra- 
dition, by demanding that a different pulse sequence be used 
for every given T: given a feasible pulse interval To, the ap- 
propriate CDD order n is found from the condition T = 4 n To. 
Unlike in PDD, the pulse order in CDD level n + 1 is com- 
pletely different from that in level n. 

Next we present the results of the first-ever CDD exper- 
iment. 13 The results demonstrate that, as predicted) 6 ! 9 ! 20 ^ 
CDD outperforms PDD in preserving quantum memory, and 
CDD becomes increasingly effective at higher concatenation 
levels. 



III. EXPERIMENTAL RESULTS ON QUANTUM 
MEMORY USING CDD 

Presented in Figures [T][3] are the results of quantum mem- 
ory NMR spectroscopy experiments. The experiments were 
performed at TU Dortmund with a homebuilt 300 MHz solid- 
state NMR spectrometer using a 7.05 T Oxford magnet and a 
homebuilt double resonance (HC) probe. Before the experi- 
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merits, we first used suitable tune-up pulse sequence d 30 " 32 ! to 
minimize pulse-length errors and phase transients for the 13 C 
channel by adjusting the probe tuning and RF amplitudes. For 
the tune-up, we used a liquid sample of 13 C-labeled Methanol 
(99% enriched). 

To measure the effect of the decoupling sequences, we used 
a sample of powdered Adamantane. In this plastic crystal, the 
nearly spherical molecules tumble rapidly and isotropically in 
the solid phase. The motion averages all intramolecular dipo- 
lar couplings to zero, but does not eliminate intermolecular 
couplings. As a result of this averaging process, there is only 
one coupling between every pair of molecules, thereby reduc- 
ing the Adamantane molecule to a point dipole source contain- 
ing 16 proton spins or 16 proton spins and one 13 C-nuclear 
spin in any C position (11% of the molecules). All remain- 
ing interactions (except for the Zeeman term) are insignificant. 
The sensitivity of 13 C NMR signals was enhanced by the stan- 
dard cross-polarization experiment, which provides polariza- 
tion transfer from the abundant 1 H spins. The enhancement 
factor of the 13 C signal was about 3.2. After the transfer, the 
carbon spin magnetization was stored for 2 ms as longitudi- 
nal magnetization. Another tt 12 pulse created transverse 13 C 
magnetization as the initial state for the dynamic decoupling 
experiments. 

In the experiments, the width S of the tt pulses for 13 C was 
10.52 [is. Pairs of tt pulses with the same phase and no inter- 
vening delay were omitted. For pairs of tt pulses with dif- 
ferent phases, we inserted a minimum delay f a = 376 ns 
to change the phase of the rf-pulse. As an example, XZf 
(where / denotes free evolution, without pulses) was imple- 
mented as Xf a Zf. The minimum cycle time for the PDD 
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- 5), for the even-order CDD se- 
2 f a , and for the odd-order CDD 
4(r^ 2fc + S), where To is the interval be- 
tween consecutive tt -pulses. The remaining signal was de- 
tected by measuring the free induction decay of the 13 C spins 
while decoupling the protons after the last echo pulse and a 
dead- time of = 8.5/iS. To obtain an accurate signal am- 
plitude, we Fourier-transformed the free induction decays and 
integrated the signals over the relevant frequency range. 

Figure [l] shows the spin echo signal for the 13 C nuclear spin 
qubits. In each experiment, we applied cyclic sequences of tt- 
pulses corresponding to the various CDD orders and measured 
the remaining signal after the end of the sequence. We plot the 
absolute value of the transverse magnetization, corresponding 
to 



(a x - i(jy) \ = J Tr[a x p(T t 



r Pk\}2 



■Tr[a y p(T t 



Pk\]2 



(7) 



where the density matrix p{r^ k ) represents the state of the 
system at the end of DD pulse sequence (either PDD or CDD, 
where k now represents the number of cycle repetitions (for 
PDD) or number of concatenations (for CDD)). Figure [2] 
shows the decay of the signal as the number of cycles and 
thus the total evolution time increases. In these experiments, 
the delay between pulses was set to To = 15/is. By fit- 
ting the experimental data to an exponential decay function 
S = S§e~ t l T ' 2 , we extracted the decoherence rates I/T2. Fig- 



ure [3] shows the decoherence rates extracted from a series of 
experiments conducted with different pulse intervals To . 

These data show that: (i) higher levels of CDD suppress 
errors more efficiently than lower levels, (ii) at high concate- 
nation levels the benefits of CDD increase as To decreases, 
(iii) CDD always and substantially outperforms PDD, an im- 
provement that becomes more pronounced as To decreases. 
Indeed, comparing the results for CDDi = PDD and CDD3 
at To = 15/is in Fig. [3] we observe an improvement of almost 
an order of magnitude in the decoherence rate. 

We note the decoherence here is not engineered, but in- 
stead caused by the actual dipole-dipole interactions between 
the 13 C qubit and the 1 H spins inherent in the Adamantane 
molecule. We further note that the specific nature of the in- 
teraction, be it dipole-dipole, Heisenberg, or something else 
altogether, is inconsequential to the success of CDD. What 
matters is the strength of this interaction relative to the DD 
pulse interval To, as will be made more explicit in the theoret- 
ical analysis, to which we now turn. 
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FIG. 2: Signal decays of 
CDD 2 and CDD 3 with t = 
to an exponential decay. 



3 C-Spin for the free evolution, CDDi, 
15/xs. The solid lines denote the fitting 



IV. THEORETICAL ANALYSIS 

We now present a theoretical interpretation of some of 
the results shown in Figures [T][3] To allow for a generic 
analysis of the DD response to different bath environments 
and system-bath couplings, we characterize the leading-order 
DD behavior for simplicity in terms of the following rele- 
vant parameters, which capture the strength or overall rate 
of the internal bath and system-bath dynamics, respectively: 
j3 = \\H B \\ and J = \\H S b\\^ The total Hamiltonian is 
H = H s + H B + HsB-If J > A then the system-bath cou- 
pling is a dominant source of error, and remains relatively sta- 
ble since the internal bath dynamics are comparatively slow. 
In this case, DD should produce significant fidelity gains as 
it removes the dominant error source. On the other hand, if 
J < f3, then the system coupling to the environment induces 
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FIG. 3: Decoherence rate vs. pulse interval for the system shown in 
Figure [1] for CDDi = PDDi (red, circles), CDD 2 (blue, squares), 
and CDD 3 (black, triangles). For PDD with non-ideal pulses, errors 
due to finite pulse-width accumulate and overwhelm the improve- 
ment due to the smaller pulse interval. CDD, on the other hand, 
compensates for finite pulse width errors, and hence improvement is 
seen with smaller pulse intervals, as expected. 

relatively slow dynamics, while the environment itself has fast 
internal dynamics. In this case, minimizing the system-bath 
coupling will have less of an effect on the overall dynamics, 
so it may be considered a worst case scenario when assessing 
DD performance. 

Let us now compare the pure system state p° s obtained in 
the error- free setting, where Hsb = (the desired state), to 
the system state ps resulting from the presence of a non-zero 
Hsb, but subject to CDD level n (the actual state). The appro- 
priate distance measure, which we would like to minimize, is 
the trace-norm distance D, defined in Appendix [A] For a total 
evolution time T = Ntq using TV = 4 n zero- width pulses, the 
distance D between the desired and actual state is bounded, in 
the "pessimistic" regime J < f3, as! 20 * 22 * 34 ! 

D[ps(T),p° s (T)}<2JT e n , (8) 

where e = 4/3ro2 n . Here e plays the role of a "threshold 
parameter", i.e., rapid convergence of the distance to zero in 
the large n limit is guaranteed if e < 1. This situation is 
reminiscent of the error threshold in quantum fault tolerance 
theory, 35 but requires us to decrease To while the concatena- 
tion level rises, in order to mitigate the 2 n factor. In general 
this is undesirable since the pulse interval cannot be decreased 
indefinitely. Indeed, in each of our experiments the pulse in- 
terval was fixed and the total time was variable. The interval 
length was varied from experiment to experiment, and the to- 
tal time was adjusted accordingly, depending on the number 
of pulses used. Therefore we next analyze the above bound on 
D under the assumption of fixed pulse interval To. We shall 
show that, perhaps counterintuitively, there is a regime where 
increasing the concatenation level, and hence also the total 
evolution time, actually leads to improved fidelity. Such an 
improvement can, however, not continue indefinitely, and we 
shall find that there is an optimal concatenation level beyond 



which fidelity starts to drop. 

If 5* is the desired maximum probability of distinguishing 
Ps(T) from p° s (T), we can set the bound ^ on D equal to 
(5*. This yields a quadratic equation in the number of concate- 
nation levels n, whose physically meaningful branch is: 

n= - vTl + log 4 (/3r )] 2 - 21og 4 (Jr /(5*) - 1 

-[l + log 4 08T )]. (9) 

(Rounding to the nearest integer is implied.) This quantifies 
the "simulation overhead" in CDD (the analog of the poly- 
logarithmic simulation overhead in fault tolerance theory^, 
and means that using CDD, in the limit of zero- width pulses, 
we can obtain a given simulation accuracy target with a 
number of concatenation levels that scales only logarithmi- 
cally in /3to, Jto, and 5*. We can also minimize the bound 
^ on D for n at fixed Jto and /3ro . The result is 

n opt = [log 4 - 1J, (10) 

which is the optimal number of concatenation levels, in agree- 
ment with more involved derivations PS At this optimal 
value, which increases with decreasing /3tq, the bound on 
D scales as e~ 2 0°S4 P T ®) and also becomes very small in 
the small /3ro limit. Since the internal bath dynamics are 
the dominant error source here, it is not surprising that the 
optimal concatenation level only depends on /3tq. A simi- 
lar expression for optimal concatenation level occurs i n the 
"optimistic" J > j3 regime, depending on both J and /3.I2EI 
The key point is that even in the simplest case: ideal, single- 
qubit quantum memory, there exists an optimal concatenation 
level depending on the Jto and /3to parameters, after which 
point increasing concatenation no longer improves system fi- 
delity. Intuitively, this optimal concatenation level represents 
the turning point between the competing processes of im- 
proved CDD error suppression and the exponentially increas- 
ing time for error-prone evolution that occurs with each level 
of concatenation. The precise point at which DD error sup- 
pression is overcome by the accumulated error depends on 
the relative strengths of the error-inducing Hamiltonians, de- 
termined by Jto and /3tq. For comparison, let us also con- 
sider the case of PDD with N pulses. In this case the bound 
D[p s (T),p° s (T)] < 2NJr f3r applies under the same as- 
sumptions as those required for the CDD bound |8| on D. 
The evidently much more rapid convergence of CDD, as long 
as pulse intervals shorter than l/(2 n+2 /3) are attainable (the 
e < 1 condition), is due to its recursive structure, which allows 
CDD to cancel successive orders of the s ystem -bath propaga- 
tor in time-dependent perturbation theorypQEH 

Returning now to our experimental CDD results (Fig. [3}, 
which are in the "pessimistic" regime J < /3, 36 we see that 
for CDD3 they agree with the theoretical prediction that un- 
der concatenation, lowering To improves the fidelity. 37 More- 
over, they are consistent with the e < 1 regime, where CDD 
outperforms PDD. However, as can be seen from Figure [T] 
while CDD consistently outperforms PDD, and both are bet- 
ter than free (unprotected) evolution, the performance of both 
PDD and CDD deteriorates with increasing total experiment 
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time T. Thus, our experiments have n opt = 1 in the sense 
of the optimal concatenation level of Eq. ( [TO] ). In other words, 
our experiments did not reach the "sweet spot" where a higher 
level of concatenation (and hence total evolution time) actu- 
ally leads to increased fidelity, which would lead to n opt > 1. 
The closeness of e to 1 and the finiteness of the pulse width 
play a role in this regard. 

V. SUMMARY AND CONCLUSIONS 

Early DD schemes were designed to remove unwanted 
system-bath interactions to a given, low order in time- 
dependent perturbation theory^R CDD is the first explicit 
scheme proposed to remove such interactions to an arbi- 
trary order, via a recursive construction in which each suc- 
cessive level removes another order in time-dependent per- 
turbation theory 6 . This work reports the first experimental 
demonstration of CDD, using NMR. As predicted, CDD out- 
performs PDD in preserving the fidelity of quantum states. 
The experimental system studied here — the 13 C spin qubit 
of Adamantane — undergoes single-axis decoherence, so that 
one might suspect that the full CDD sequence we used is 
"overkill" for this system, because this sequence was designed 
to suppress full three-axis decoherence. However, pulse er- 
rors and magnetic field inhomogeneity result in second order 
off-axis contributions to decoherence, which is why DD se- 
quences using just a single pulse type (e.g., UDD), cannot be 
expected to perform well on our system. Nevertheless, it is 



important to identify an experimentally feasible NMR system 
where dephasing and longitudinal relaxation are of compa- 
rable magnitude, and to study the performance of CDD rela- 
tive to other methods capable of suppressing full decoherence, 
such as QDD. 
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Appendix A: Distance and fidelity 

A distance measure between quantum states p and a 
(density matrices) is provided by the "trace-norm distance": 
D(p,v) = !||p-<t||i, where ||A||i = W^U = Si (A), 
and where Si(A) are the eigenvalues of VA^A. The trace 
norm distance is the maximum probability of distinguishing p 
from a, vanishes if and only if p = a, and is related to the 
fidelity F(p,a) = HVpV^IIi via 1 - D(p,a) < F(p,a) < 
yl — D(p, a) 2 , so that knowing one bounds the other. In the 
case of comparing an output state a to a desired pure state 
p — (ip\, this reduces to F = \J a \ifj), and perfect 
fidelity (F = 1) results if and only if p = a. 
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lation: a spin system consisting of one 13 C and six 1 H with these 
parameters produced the best fit to the experimental free induction 



decay. 

We are taking some liberty here with the term "fidelity": our the- 
oretical analysis indeed provides a fidelity bound, while the ex- 
perimental results report different, but related quantities, namely 
NMR signal and decoherence rate. Moreover, our experiment is 
restricted to testing the performance of CDD for a single input 
state, whereas the theoretical distance and fidelity measures we 
have discussed apply for general input states. 



